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Abstract 

Broad line radio galaxies (BLRGs) arc a rare type of radio-loud AGN, in which the broad optical 
permitted emission lines have been detected in addition to the extended jet emission. Here we report on 
deep (40ksec x4) observations of the bright BLRG 3C 120 using Suzaku. The observations were spaced 
a week apart, and sample a range of continuum fluxes. An excellent broadband spectrum was obtained 
over two decades of frequency (0.6 to 50 keV) within each 40 ksec exposure. We clearly resolved the iron 
K emission line complex, finding that it consists of a narrow K Q core (er ~ 110 cV or an EW of 60 eV), 
a 6.9 keV line, and an underlying broad iron line. Our confirmation of the broad line contrasts with the 
XMM-Newton observation in 2003, where the broad line was not required. The most natural interpretation 
of the broad line is iron K line emission from a face-on accretion disk which is truncated at ~ 10 r g . Above 
10 keV, a relatively weak Compton hump was detected (reflection fraction of R ~ 0.6), superposed on 
the primary X-ray continuum of T ~ 1.75. Thanks to the good photon statistics and low background of 
the Suzaku data, we clearly confirm the spectral evolution of 3C 120, whereby the variability amplitude 
decreases with increasing energy More strikingly, we discovered that the variability is caused by a steep 
power-law component of T ~ 2.7, possibly related to the non-thermal jet emission. We discuss our findings 
in the context of similarities and differences between radio-loud/quiet objects. 
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1. Introduction 

One of the most important issues in the study of active 
galactic nuclei (AGN) is why well-collimated, powerful, 
relativistic radio jets exist only in 10 % of the AGN class, 
i.e., in so-called radio- loud objects (e.g., Urry & Padovani 
1995). 1 All AGNs are thought to be powered by accretion 
of matter onto a supcrmassive black hole, presumably via 
an equatorial accretion disk. Recent VLBI observations 
of a nearby active galaxy M87 confirmed that the jet is 
already launched within ^60 r g (where r g = GM/c 2 is 
the gravitational radius), with a strong collimation occur- 

1 More accurately, a well-known dichotomy of the radio-loud and 
radio-quiet objects is based on the old VLA study of PG quasars 
(Kellermann et al. 1989), and recent deep VLA FIRST survey 
shows no such sign of bimodality in radio loudness (White et al. 
2000). 



ring within ^200 r g of the central black hole (Junor et 
al. 1999). These results are consistent with the hypoth- 
esis that jets are formed by an accretion disk, which is 
threaded by a magnetic field. Therefore the observational 
properties of the accretion disk and corona are essential 
ingredients to jet formation (e.g., Livio 1999 and refer- 
ences therein). 

The profile of the iron K Q (6.4 keV) line can be used 
to probe the structure of the accretion disk, because it is 
thought to result from fluorescence of the dense gas in the 
geometrically thin and optically thick regions of the inner 
accretion disk (~ 10 r g ). The most famous example is the 
spectrum of the Seyfert 1 (Sy-1) galaxy MCG-6-30-15, 
which shows a relativistically broadened Fe K Q emission 
line, first detected by ASCA (Tanaka et al. 1995; Iwasawa 
et al. 1996; 1999). This finding was confirmed by Chandra 
and XMM-Newton (e.g., Wilms et al. 2001; Fabian et al. 
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2002), and most recently by Suzaku (Minuitti et al. 2006). 
Similar broad relativistic iron line profiles have been de- 
tected in several other typc-1 AGN (e.g., Nandra et al. 
1999, Iwasawa et al. 2004, Turner et al. 2005 for NGC 
3516; Page et al. 2001 for Mrk 766), although they are 
perhaps somewhat less common than anticipated from the 
ASCA era (e.g., Nandra et al. 1997). In contrast, the pres- 
ence of a narrower 6.4 keV line from more distant matter 
(e.g., from the outer disk, broad/narrow line regions, or 
torus) is common in many type-1 AGNs (e.g., Yaqoob 
& Padmanabhan 2004). In this context, studies of the 
iron line profile in radio-loud AGN provides important 
clues to the disk-jet connection, particularly by compar- 
ison with Sy-ls. Moreover, in a standard picture of two- 
phase disk-corona model proposed by Haardt & Maraschi 
(1991), the corona temperature is related to the high en- 
ergy cutoff typically observed in the hard X-ray band (see 
also Poutanen & Svennson 1996). Thus observations of 
the Compton reflection hump and its cutoff above 10 keV 
can be used to determine the geometry and temperature 
of the disk-corona system postulated to exist in both types 
of AGNs (e.g., Wozniak et al. 1998). 

3C 120 (z = 0.033) is the brightest broad line radio 
galaxy (BLRG), exhibiting characteristics intermediate 
between those of FR-I radio galaxies and BL Lacs. It has a 
one-sided superluminal jet on 100 kpc scales (Walker et al. 
1987; Harris et al. 2004), and superluminal motion (with 
an apparent velocity /? app = 8.1) has been observed for 
the jet component (Zcnsus 1989). This provides an upper 
limit to the inclination angle of the jet to the line of sight 
of 14 deg (Eracleous & Halpern 1988). Interestingly, the 
optical spectrum of 3C 120 is not LINER-like, as is often 
seen in FR-I radio galaxies (e.g., Baum et al. 1995), but 
rather typical of Sy-ls. It resides in an optically peculiar 
galaxy that shows only some indication of spiral structure 
(Moles et al. 1988). From reverberation mapping, the cen- 
tral black hole mass is remarkably well constrained: M = 
5.5±l \ x 10 7 M Q (Peterson et al. 2004; see also Wandel, 
Peterson & Malkan 1999). Very recently, Marscher et al. 
(2002) found that dips in the X-ray emission of 3C 120 are 
followed by ejections of bright superluminal knots in the 
radio jet, which clearly indicates an important connection 
between the jet and the accretion disk. 

In X-rays, 3C 120 has been known to be a bright (~ 
5x 10 -11 erg cm~ 2 s _1 at 2—10 keV), variable source with 
a canonical power-law spectral shape that softened as the 
source brightened (e.g., Maraschi et al. 1991 and refer- 
ence therein). A broad iron K a line was first detected by 
ASCA in 1994, with its width a = 0.8 keV and EW (equiv- 
alent width) ~400 eV (Grandi et al. 1997; Reynolds 1997; 
Sambruna et al. 1999). The follow-up observations by 
RXTE (Eracleous, Sambruna & Mushotzky 2000; Gliozzi, 
Sambruna & Eracleous 2003) and BeppoSAX (Zdziarski 
& Grandi 2001) also detected a broad iron line, but with 
a much smaller EW of ~100 eV. Furthermore, these ob- 
servations have confirmed a presence of a weak Compton 
hump in 3C 120, with a reflection normalization of f2/2-7r 
~ 0.4. It was argued that both the weak line and rela- 
tive weakness of the Compton hump is suggestive of an 



optically thick accretion disk which transitions to a hot, 
optically thin flow (Eracleous, Sambruna & Mushotzky 
2000; Zdiarski k Grandi 2001). Without a doubt, 3C 120 
is a key source to construct a unified view of radio-loud 
and radio-quiet objects. 

Most recently, 3C 120 was observed for nearly a full or- 
bit (130 ksec) with XMM-Newton on 26-27 August 2003 
(Ballantyne et al. 2004, Ogle et al. 2005). This clearly 
confirmed the presence of the neutral Fc line emission 
(57±7 eV in EW), which was slightly broadened with a 
FWHM of a = 9000±3000 km s" 1 . Both of these papers 
argued that the line profile is rather symmetric and no 
evidence was found for relativistic broadening, or alterna- 
tively the line arises from an accretion disk radius of >75 
r g at an inclination angle of — 10 deg (where relativistic 
gravitational effects are almost negligible) . They also con- 
firmed a weaker emission line at 6.9 keV with EW = 20±7 
eV (initially suggested by Chandra HETGS; Yaqoob & 
Padmanabhan 2004) , which can be interpreted as a blend 
of Fe K/3 and H-like (or He-like) iron lines. Despite signif- 
icant progress made by XMM-Newton, most of the results 
are not conclusive due to XMM-Newton's limited energy 
range (no coverage above 10 keV), and the relatively high 
background of XMM-Newton above 5 keV. At present, 
various models can fit these iron lines equally well, though 
they assume quite different geometries and/or ionization 
states (Ballantyne et al. 2004, Ogle et al. 2005). It has 
been argued that such degeneracies may be resolved by 
Suzaku, due to its unprecedented sensitivity between 0.3 
and -300 keV. 

In this paper, we present a detailed analysis of 160 
ksec worth of data on 3C 120, observed with Suzaku in 
February and March 2006 as a part of the SWG (science 
working group) program. As our aim was to monitor the 
source in various states of source activity, we divided this 
total exposure into four exposures of 40 ksec each, with 
one pointing per week (to be roughly equal to the vari- 
ability timescale of this source; e.g., Marshall et al. 2003; 
Gliozzi et al. 2003). Thanks to the excellent energy reso- 
lution and sensitivity of the XIS and HXD /PIN onboard 
Suzaku, we successfully obtained one of the highest quality 
data on this radio galaxy ever reported, between 0.6 and 
50 keV. The paper is organized as follows. The observa- 
tion and analysis methods are described in §2. We present 
an overview of 3C 120's variability during the Suzaku 
observations and temporal studies of light curves in §3. 
Detailed spectral studies are presented in §4; in particu- 
lar we focus on (1) multiband spectral features, (2) the 
nature of the iron K line complex, and (3) the difference 
spectrum between high and low states. Based on these 
new findings, we discuss the nature of the variability and 
spectral evolution of 3C 120 in §5. Finally, we present 
our conclusions in §6. Uncertainties of background mod- 
els on the HXD/PIN light curve will be discussed further 
in detail in the appendix. 



Suzaku Observations of 3C 120 

Table 1. Suzaku observation log of 3C 120. 
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obsJD start (UT) stop (UT) XIS/HXD cxp. (ksec) 

3C 120 #1 2006 Feb 09 03:20 2006 Feb 10 05:50 33.3/29.3 
3C 120 #2 2006 Feb 16 13:08 2006 Feb 17 14:06 36.4/32.6 
3C 120 #3 2006 Feb 23 20:02 2006 Feb 24 18:00 37.0/36.2 
3C 120 #4 2006 Mar 02 22:29 2006 Mar 03 20:39 37.9/37.4 



2. Observation and Data Reduction 

The broad line radio galaxy 3C 120 was observed with 
Suzaku (Mitsuda et al. 2006) four times in February and 
March 2006 with a total (requested) duration of 160 ksec. 
Table 1 summarizes the start time and the end time, 
and the exposures of each observation. Suzaku carries 
four sets of X-ray telescopes (Serlemitsos ct al. 2006) 
each with a focal-plane X-ray CCD camera (XIS, X-ray 
Imaging Spectrometer; Koyama et al. 2006) that is sen- 
sitive in the energy range of 0.3—12 keV, together with 
a non- imaging Hard X-ray Detector (HXD, Takahashi et 
al. 2006; Kokubun et al. 2006), which covers the 10-600 
keV energy band with Si PIN photo-diodes and GSO scin- 
tillation detectors. Three of the XIS (XIS 0, 2, 3) have 
front-illuminated (FI) CCDs, while the XIS 1 utilizes a 
back-illuminated (BI) CCD. The merit of the BI CCD 
is its improved sensitivity in the soft X-ray energy band 
below 1 keV. In all four observations (3C 120 #1 — 4), 
3C 120 was focused on the nominal center position of the 
XIS detector. 

2.1. XIS Data Reduction and Analysis 

For the XIS, we analyzed version 0.7 of the screened 
data (Fujimoto et al. 2006) provided by the Suzaku team. 
The screening of the version 0.7 data are based on the 
following criteria; (1) only ASCA-grade 0,2,3,4,6 events 
(Yamaguchi et al. 2006) are accumulated, while hot and 
flickcrling pixels were removed from the XIS image using 
the CLEANSIS script (2) the time interval after the pas- 
sage of South Atlantic Anomaly (T_SAA) is larger than 
436 sec, (3) the object is at least 5 deg and 20 deg above 
the rim of the Earth (ELV) during night and day, respec- 
tively. In addition, we also select the data with the cutoff 
rigidity (COR) larger than 6 GV. The XIS events were 
extracted from a circular region with a radius of 4.3' cen- 
tered on the source peak, whereas the background was ac- 
cumulated in an annulus with its inner and outer radii of 
4.9' and 6.6', respectively. For the 3C 120 #3 observation 
(see Table 1), we discard the first 600 sec data after the 
maneuver because the pointing fluctuated from its center 
(i.e., XIS nominal position) by ~ 2'. 

The XIS spectra were corrected for the hydrocarbon 
(C24H38O4) contamination on the optical blocking fil- 
ter, by including an extra absorption column due to 
Carbon and Oxygen in all the spectral fits. The col- 
umn densities for each detector were calculated based 
on the date of the observation using an empirical rela- 
tion released by the Suzaku instrumental team. We as- 
sumed different Carbon and Oxygen column densities (Ac 



and No) for different sensors and observations; Ac = 
(2.38-2.56)xl0 18 , (3.20-3.50)x 10 18 , (3.65-3.80) xlO 18 , 
and (5.44-5.65) xlO 18 atoms cm" 2 for XIS 0, 1, 2, 3 
respectively, with the ratio of C/O column densities 
(Nq/No) set to 6. This additional soft X-ray absorption 
due to the hydrocarbon contamination was included as 
a fixed spectral component using the VARABS absorption 
model in all the spectral fits. 

Although 3C 120 is known to have an large extended X- 
ray jet that is ~20 arcsec apart from the nucleus (Harris 
et al. 2004), and could not be resolved with Suzaku, the 
X-ray contribution from this large-scale (10 kpc scale) jet 
is less than 1 % of the nucleus emission in the 0.5—10 keV 
bandpass (Fj ot ~ 2.3xl0 -13 erg cm -2 s" 1 ). We therefore 
use the latest version of both the response matrix and the 
point spread function (PSF) released by the Suzaku team 
for the point-like sources, ae_xi[0,l,2,3]_20060213.rmf and 
ae_xi[0,2,3]_xisnom6_20060415.arf. Since the nuclear emis- 
sion is very bright from 0.5—10 keV, we binned the XIS 
spectra to a minimum of 200 counts per bin to enable the 
use of the \ 2 minimization statistic. 

Another note on the XIS data analysis is the accuracy of 
the energy scale (reported to be less than 0.2 %; Koyama 
et al. 2006) and energy resolution. These calibrations are 
very important because one of our ultimate goals is to ob- 
tain the precise measurement of the iron line complex, as 
we will see in § 4.3. For this purpose, the 55 Fe calibration 
source located on the corners of the XIS chips were used as 
an accurate calibrator of the instrumental response dur- 
ing 3C 120 observations. 2 The 55 Fe source produces 
a characteristic X-ray line from Mn K Q at 5.895 keV (a 
combination of K Q i and K Q 2 at 5.899 keV and 5.888 keV 
respectively with a branching ratio of 2:1). Figure 1 shows 
the results of Gaussian fitting for each of the XIS chips 
during 3C 120 observation #1 — 4. 

From measuring the 5.895 keV lines in the calibration 
source, we find that the line energy is shifted slightly red- 
wards by 5.5±0.6 eV, while there is some residual width 
(after the instrumental response function has been ac- 
counted for) in the calibration lines of cr ca i = 45±1 eV. 
A further check of the line width was made by using Sgr 
C (molecular cloud) data which was taken between 3C 120 
#2 and #3. As reported in Ueno et al. (2006), the 1 a 
width of the 6.4 keV line in the Sgr C observation was 
39±5 eV, which is consistent with the 3C 120 calibra- 

2 Note the 55 Fe calibration source is on the corners of the XIS 
chips, while 3C120 is focused on the center of the CCD. From 
detailed analysis of galaxy clusters and Sgr C, the XIS team 
confirmed this position difference makes only negligible effects 
in both the line width and spatial gain variations. 
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Fig. 1. Energy scale and residual width of the 55 Fe calibra- 
tion source on each of the XIS sensors, during each of the four 
observations. 

tion lines within the statistical error. Since the intrinsic 
line width is expected to be negligible for this molecular 
cloud, it must be mostly instrumental. Therefore in the 
following, the intrinsic width of the iron K Q line is sim- 



ply evaluated as a- n 
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obs 



°7ai> where cr b s is the 



measured width of the iron line discussed below. 
2.2. HXD Data Reduction and Analysis 

The source spectrum and the light curves were ex- 
tracted from the cleaned HXD/PIN event files (version 
1.2) provided by the Suzaku team. The HXD/PIN data 
are processed with the screening criteria basically the 
same as those for the XIS, except that ELV > 5 deg 
through night and day, T_SAA > 500 sec, and COR > 
8 GV. The HXD/PIN instrumental background spectra 
were generated from a time dependent model provided 
by the HXD instrument team for each observation. The 
model utilized the count rate of upper discriminators and 
COR values as the measure of cosmic-ray flux that pass 
through the Si PIN diode; background spectra based on 
a database of non X-ray background observations made 
by the PIN diode to date are provided (See Kokubun et 
al. 2006; Fukazawa et al. 2006 for more details). At the 
time of writing, two different background models, A and 
B, are under investigation, but the difference between the 
two models is rather small for our 3C 120 HXD/PIN data 
analysis. According to careful examinations using simula- 
tions on various observational datasets, the uncertainty of 
the background models for the PIN detector is expected 
to be ~ 5% for both model- A and B (see Appendix). 

Both the source and background spectra were made 
with identical GTIs (good time intervals) and the ex- 
posure was corrected for detector deadtime of ~ 5% 
(ranges in 3.5—8.5 % for 3C 120 observations; sec 
Appendix). To minimize uncertainties due to limited 
photon statistics on the background models, background 
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Fig. 2. The combined HXD/PIN spectra for the Suzaku ob- 
servation of 3C 120 #1-4 over the whole HXD/PIN en- 
ergy bandpass (10—60 keV). Black shows the total spectrum 
(source plus background) , blue shows sum of non- X-ray back- 
ground and the CXB, and red for the net source spectrum. 
Note the source is detected in each data bin at more than 10 
% of background level over the 12—50 keV range used in this 
paper. ) 

spectra were generated with 10 times the actual back- 
ground count rate but increasing the effective exposure 
time of background by a factor of 10. The HXD/PIN re- 
sponse file dated 2006-08-14 for the XIS nominal position 
(ae_hxd_pinxinom_20060814.rsp) was used for these spec- 
tral fits. Full details of the HXD instrument and perfor- 
mance are given in Takahashi et al. (2006) and Kokubun 
et al. (2006). 

The time averaged (obs #1-4 combined) HXD/PIN 
spectrum thus obtained is shown in Figure 2, plotted over 
the energy range of 10 to 60 keV. HXD/PIN data be- 
low 12 keV have been ignored to avoid noise contamina- 
tion near the lower threshold of the PIN diode. Also the 
data above 50 keV are discarded, as a detailed study of 
noise and background systematics is on-going above this 
energy. Figure 2 shows the total (3C 120 + observed 
background) spectrum, where the background includes 
both the instrumental (non X-ray) background model-A 
and the contribution from the cosmic X-ray background 
(CXB) (Gruber et al. 1999). Here the form of the CXB 
was taken as 9.0x 10~ 9 (E/3keV)-°- 29 cxp(-E/40kcV) erg 
cm -2 s" 1 sr" 1 keV -1 and the observed spectrum was sim- 
ulated assuming the PIN detector response to isotropic 
diffuse emission. When normalized to the field of view of 
the HXD/PIN instrument the effective flux of the CXB 
component is expected to be 9.0xl0 -12 erg cm~ 2 s _1 in 
the 12-50 keV band, which is about -14 % of the 3C 120 
flux in the same energy bandpass. 

3C 120 is known to have a 50-100 keV flux of ~ 
5xl0 -11 erg cm 2 s -1 (e.g., Wozniak et al. 1998; Zdziarski 
& Grandi 2001) and therefore ultimately can be detected 
by the HXD/GSO detector. However, this is only a few 
percent of the GSO detector background and the study 
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Fig. 3. The overall variability of 3C 120 during February and 
March observations by Suzaku. upper: 0.4—2 keV XIS (XIS 
0-3 summed), middle: 2-10 keV XIS (XIS 0-3 summed), 
and lower: 12—40 keV HXD/PIN. The backgrounds are sub- 
tracted. 

at this level of background systematics is still on-going by 
the HXD instrumental team. Therefore in this paper, we 
do not include the GSO data in the subsequent spectral 
fits. Above 50 keV, a reliable detection of 3C 120 cannot 
be made at the present time by the GSO, but results us- 
ing new response matrices and revised background models 
will be discussed elsewhere in the near future. 

3. Temporal Studies 

3.1. Overview of Variability 

Figure 3 shows an overview of the count rate variations 
of 3C 120 during the February and March observations 
(#1-4). The light curves of the 4 XISs and PIN detectors 
are shown separately in different energy bands; 0.4—2 keV 
(upper; XIS), 2-10 keV (middle; XIS) and 12-40 keV 
(lower; HXD/PIN). The net source count rates, averaged 
over four observations, measured from 0.4—10 keV were 
2.965±0.005 cts s~\ 3.703±0.005 cts s" 1 , 2.972±0.005 cts 
- 1 , and 2.799±0.005 cts ~\ respectively for the XIS 0, 1, 
2,3. The background is typically 3 % of the source counts 
for the FI-XIS (XIS 0,2,3) and 7 % for the BI-XIS (XIS 1). 
For the PIN detector, the net average source count rate in 
the 12-40 keV band was 0.150±0.002 cts _1 , compared 
to the PIN background (non X-ray background + CXB) 
rate of 0.449 cts _1 . 

Figure 3 clearly indicates that 3C 120 was in a relatively 
high state during the 1st 40 ksec observation (#1; the 
summed count rates of 4 XISs detectors was 15.94±0.01 
counts s" 1 ), then its count rate dropped by ~20% in the 
2nd observation, and finally reached a minimum in the 
4th observation (#4; 12.02±0.01 counts s^ 1 ). Here the 
effect of time dependent degradation of the XIS efficiency 
due to C and O contamination (Koyama et al. 2006; see 
also §2.1) is corrected between 3C 120#1 and #4, though 
this makes only negligible effect even below 2 keV within 
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Fig. 4. (a) upper: Correlation of XIS count rates between 
0.4—2 keV and 2—10 keV. (b) lower: Changes of hardness 
ratio between 0.4—2 keV and 2—10 keV. The hardness is de- 
fined as 2—10 keV count rates divided by 0.4—2 keV count 
rates. 

a month scale. Interestingly, a large time variation of the 
0.4—2 keV count rate in the 1st observation was not clearly 
detected above 2 keV, suggesting a relative lack of vari- 
ability at higher energies, as has already been noticed for 
this particular source by various authors (e.g., Maraschi 
et al. 1991 and Zdziarski & Grandi 2001). The count rate 
variations of the HXD/PIN detector is less clear due to 
limited photon statistics, but some level of variability ex- 
ists, as we will see in the next section. 

3.2. Hardness Ratio 

Figure 4 (a) compares the XIS count rate correlation 
between the soft X-ray (0.4—2.0 keV) and hard X-ray en- 
ergy bands (2—10 keV). Generally the correlation is tight, 
such that hard X-ray flux increases when the soft X-ray 
flux increases. However, some scatter in the correlation 
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Fig. 5. (a) upper: Correlation of count rates between 2—10 
keV (XISs) and 12-40 keV (HXD/PIN). (b) lower: Changes 
of hardness ratio between 2—10 keV and 12—40 keV. The 
hardness is defined as the 12—40 keV count rate divided by 
the 2—10 keV count rate. 

suggests some variation of the spectrum, even if the source 
is in similar flux states. Moreover, the amplitude of flux 
variations is approximately 40 % in 0.4—2 keV, but only 
~ 25 % in 2-10 keV. Figure 4 (b) shows the soft X-ray 
count rates (0.4—2 keV) versus hardness ratio, where the 
hardness is defined as 2—10 keV count rates divided by 
0.4—2 keV count rates. This clearly suggests a spectral 
evolution such that the spectrum becomes softer when 
the source becomes brighter, but again the correlation is 
rather loose, especially when the source is in lower states 
of activity (obs #2 — 4). Also, it appears that the hardness 
diagram cannot be fit linearly against 0.4—2 keV source 
counts. This may suggest more than just one spectral 
component is responsible for the 3C 120 X-ray emission. 

Similarly Figure 5 (a) compares the 2—10 keV XIS 
and 12-40 keV HXD/PIN count rates during the 3C 120 
#1 — 4 observations. Again, one can find a weak COrrela- 
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Fig. 6. Energy dependence of variability of 3C120 #1 — 4. 
The variability parameter, excess variance, was calculated for 
the total exposure in 7 (6 for 4 XISs + 1 for HXD/PIN) energy 
bands. The HXD/PIN data where the background is high 
(photon counts of > 0.7 counts s — 1 ) are not used to minimize 
uncertainty due to background subtraction (sec Appendix). 



tion between source variability in these two energy bands. 
In this particular case, however, the hardness seems to 
stay almost constant within statistical errors as shown in 
Figure 4 (b). The apparent small changes in hardness ra- 
tio may suggest a relative lack of spectral evolution above 
2 keV, though hardness changes of ~20 % level (as ob- 
served in Figure 4 (a)) are difficult to detect above 10 
keV due to limited photon statistics. 

3.3. Excess Variance 

For more detailed temporal studies, obs #1 — 4 light 
curves are further divided into 6 energy bands (0.4—1 keV, 
1-2 keV, 2-3 keV, 3-5 keV, 5-7 keV, 7-10 keV) for the 
XISs and 1 energy bands (12-40 keV) for the HXD/PIN, 
respectively. To estimate the amplitude of variability in 
a systematic way, "excess variance" (e.g., Zhang et al. 
2002 and reference therein) was calculated for light curves 
derived in different energy bands. Excess variance (of ms ), 
is a net variance, which is defined as the difference between 
total variance (of ot ) and noise variance (c^ oisc ), that is 
caused by statistical errors; 



1 



iV 



Nx 2 f-f 



[(Xi - x) 



X 



(1) 



where Xi is the z-th bin in the light curve and x is the mean 
of Xi. The error on c^ms i s estimated by srj/(x 2 ^/N), 



where sn is the variance of quantity 
given by (Turner et al. 1999) 
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and 



(2) 



The fractional variability parameter F va r used in this pa- 



per is the square root of excess variance: F var = (<7 lr 
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Figure 6 shows F var measured in this way, using the 
overall light curves combined in obs #1 — 4. Again, we 
carefully consider the effect of time dependent degrada- 
tion of the XIS efficiency, but it only accounts for ^2 
% decrease of the XIS count rates even in 0.4—1.0 keV 
band. Clearly, the variability is larger in the lower energy 
bands; the largest variability was observed in the 0.4— 
IkeV band with F var = 14.1±1.1 %, and it gradually de- 
creases with increasing energy, and reaches 4.8±0.5 % in 
7—10 keV. The variability amplitude above 10 keV cannot 
be well constrained (F var = 5.5±3.4 %), but is consistent 
with those of 7—10 keV variability observed with the XIS. 
Note this range of variability is consistent with that re- 
ported in Markowitz & Edelson (2004), who find F var = 
6 % on 6-day time scales and 8 % on 36-day time scale in 
2—12 keV. We also calculated the excess variance of the 
light curves accumulated within each of the 40 kscc obser- 
vations (obs #1, 2, 3, 4, respectively). The amplitude is 
smaller on this shorter time scale, with F var < 3 % in all 
energy bands, and sometimes is consistent with no vari- 
ability. We therefore conclude that the variability time 
scale of 3C 120 is much larger than one day, typically t var 
~ 10 6 (sec). This relatively long timescale of variability is 
consistent with those found in literature (e.g., Marshall ct 
al. 2003; Gliozzi ct al. 2003; also see the Appendix for ap- 
parent short-term variations seen in the HXD/PIN light 
curve). 

4. Spectral Studies 

4-.1. Overview of the Broad-band Spectrum 

Before going into the detailed spectral study, we quickly 
overview the time-averaged spectrum of 3C 120 between 
0.6 and 50 keV, from the whole Suzaku observation (obs 
#1-4 combined). The XIS and HXD/PIN background 
subtracted spectra were fitted using XSPEC vll.3.2p, 
including data over the energy range 0.6—50 keV. The 
Galactic absorption toward 3C 120 is taken to be Ah = 
1.23xl0 21 cm~ 2 (Elvis, Wilkes & Lockman 1989). All 
errors are quoted at 68.3% (lcr) confidence for one inter- 
esting parameter unless otherwise stated. All the fits in 
this section are restricted in the energy range of 0.6—12 
keV (XIS 0,2.3: FI chips), 0.6-10 keV (XIS 1: BI chip), 
and 12-50 keV (HXD/PIN). 

As a first step, we fitted the 4 XIS spectra using a sim- 
ple absorbed power-law model in the 3—10 keV band ex- 
cluding 5—7 keV, in order for the continuum parameters 
not to be affected by possible iron line features, absorp- 
tion/excess features in the lower energy band, as well as 
the reflection component expected to be present above 10 
keV. Since we found that the spectral parameters of the 
4 XIS spectra are all in good agreement within the er- 
ror bars (at the few percent level), data from 3 FI-XISs 
(XIS 0, 2, 3) were co-added to maximize the signal to 
noise. Figure 7 shows the 4 XISs + HXD/PIN spectra 
with residuals to this baseline model, where the XIS data 
around the Si K edge from 1.7—1.9 keV were discarded in 
all 4 XIS chips, due to uncertain calibration at this energy 
at the time of this writing. We also did not use the data 
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Fig. 7. The broad-band (0.6-50 keV; 3 FI-XISs + PIN) 
Suzaku spectrum of 3C 120. The upper panel shows the data, 
plotted against an absorbed power-law model of photon index 
T = 1.74 and column density 1.23xl0 21 cm~ 2 , fitted over the 
4—12 keV band. The lower panel shows the data/model ra- 
tio residuals to this power-law fit. Deviations due to (1) the 
iron K-shell band, (2) the Compton reflection hump, and (3) 
excess soft emission arc clearly seen. 

below 0.6 keV, as uncertainties in C/O contamination are 
more significant (see §2.1). 3 

The best fit spectral power-law index thus determined 
was r = 1.74, but was statistically unacceptable if we ex- 
trapolate the model to lower and higher energies (x 2 /d.o.f 
= 6593/614). The residuals of Figure 7 shows that the 
spectrum exhibits features at different X-ray energies; (1) 
iron line features around 6 keV, (2) a hard X-ray bump 
above 10 keV, and (3) excess emission below 3 keV. In 
particular, the structure around the iron K emission line 
is rather more complicated than that inferred from any 
previous X-ray satellite. A zoom-in of the "iron-line com- 
plex" observed with Suzaku is presented respectively for 
the FI-XISs and BI-XIS in Figure 8. Note the asymmetric 
line profile and the presence of a red-tail below 6 keV. 

4-2. The Baseline Continuum Emission 

To model the overall X-ray spectrum between 0.6 and 
50 keV, we start from the shape of the baseline continuum 
emission. As for many radio-quiet AGNs, it has been sug- 
gested that the continuum of 3C 120 is composed of a 
direct (primary) power-law plus reflection by cold mat- 
ter surrounding the nucleus (e.g., Wozniak et al. 1998; 
Zdziarski & Grandi 2001). We assumed a primary com- 
ponent including an exponential cut-off at high energies 
of the form cx E~ r x e ~ E / E!old , where £foid is the cut- 
off energy in keV, T is the differential spectral photon 
index. Furthermore a reflection component produced by 



We note, an uncertainty in calibration affects only the "spectral 
shape" around the Si K edge. Similarly, the XIS spectra below 
0.6 keV is still uncertain due to the C/O contamination, but no 
problems arises for the temporal studies as we have presented 
in §3. 
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Table 2. Results of spectral fits to the 3—50 keV XISs + PIN time-averaged, co-added continuum spectrum of 3C 120. The 4 XISs 
data below 3keV (soft excess) and between 5—7 keV (Fe-line complex) were discarded to avoid the complexity of the models. 



Model T 


[kcV] 


R 

[£l/2ir] 


Abund 


i 

[deg] 


TP b 
^3-10keV 


-FlO-50keV C 


xVd.o.f 


PRV a 1.73t££j 


mo* 


-0.12 


1.0* 


18.2* 


30.0±0.1 


62.1±0.6 


1112/1060 



a PEXRAV model in XSPEC (agdziarz & Zdziarski 1995). We fixed the cutoff power- law energy -EWd = 100 keV and iron 
abundance Af c = 1.0. 

b 3-10 keV flux in units of 10~ 12 erg cm' 2 s -1 . 

c 10—50 keV flux in units of 10~ 12 erg cm~ 2 s _1 . Constant fraction factor between HXD and XIS was set to be Npw/Nxis = 

1.114±g;g^ (see also Figure 9). 

' Parameters fixed to these values. 
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Fig. 8. A Zoom-in of the iron line profile of 3C 120, plotted 
as a ratio against a power-law of photon index Y = 1.74. 
(upper): XIS 0,2,3 combined, (lower): XIS 1. 

Compton down-scattering of X-rays off neutral material 
was included, using the PEXRAV model within XSPEC 
(Magdziarz & Zdziarski 1995). The inclination angle i 
was fixed at 18.2 deg, which is the minimum allowed by 
the models, and also close to the limit derived from the 
superluminal motion of jet; see § 1. The solid angle f2 
subtended by the Compton reflector was allowed to vary, 
and determined by the parameter R = 0,/2-k. 

In this modeling, we fixed the Fe abundance at 
unity, following approaches previously made for this ob- 
jects in literature (e.g., Ballantync et al. 2004; Ogle et 
al. 2005). Similarly, we assumed a fixed value of -Efoid 
= 100 keV because the HXD/PIN is not sensitive above 
60 keV, and therefore cannot determine the high energy 
end of direct power-law component (in this context, see 
Wozniak et al. 1998 who confirmed a cutoff with Sfoid ~ 
100 keV at a significance of 99.95% by combined ASCA 
and OSSE analysis. Ballantyne et al. (2004) assumed a 
fixed E iold of 150 kcV to fit a combined XMM-Newton 
and RXTE/HEXTE spectrum. Also Zdziarski & Grandi 
(2000) found Sfoia ~ 120-150 keV using the Beppo-SAX 
data). 
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Reflection Fraction (R) 

Fig. 9. Contour plot showing the reflection parameter R ver- 
sus the constant normalization factor between the HXD /PIN 
and XIS, which shows that the constant factor (defined here 
as A^jjxd /pin/^xiSi where N is the normalization) is close to 
1.1, which is consistent with the cross calibration between the 
XIS and the HXD using the Crab nebula. Confidence levels 
in the figure correspond to the 68%, 90% and 99% levels of 
significance, respectively. 

As a first step, we fit the combined XIS + HXD/PIN 
data only above 3 keV, as the data below this energy is 
significantly contaminated by the soft excess emission as 
we have shown in Figure 7. Similarly, data between 5 
and 7 keV was ignored to avoid the region around the 
iron line complex. Even with this simplified model and 
constrained parameters, special care must be paid to the 
cross calibration of the XIS and the HXD/PIN, because 
it could easily affect the determination of the reflection 
parameter R. At the date of writing (2006 November), 
the relative normalization of the HXD /PIN was reported 
to be 1.13 times larger than the XIS, as determined by 
the XRT team using the spectrum of the Crab nebula 
(M. Ishida; JX-ISAS-Suzaku-mcmo 2006-40). 

Figure 9 shows the confidence contours of reflection 
parameter, R = Q/2n, against the relative normaliza- 
tion of the XISs/PIN. An apparent anti-correlation be- 
tween the two parameters indicates R is strongly af- 
fected by the normalization factor between the two in- 
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struments. Nevertheless, the mean value of normaliza- 
tion factor (1.114^ra) 

is perfectly consistent with that 
determined from the Crab nebula, and the reflection pa- 
rameter is well constrained at R — 0.47^o'i2- The re- 
sultant parameters for the baseline continuum emission 
are summarized in Table 2. The photon index of the di- 



rect component is determined to be T = 1.731q]o|. The 
observed flux of the continuum emission (direct + reflec- 
tion) is (30.0±0.1)xlCr 12 erg cm" 2 s" 1 over 3-10 keV, 
and (62.1+0.6) xlO" 12 erg cm" 2 s" 1 over 10-50 keV. 

Finally, we comment on how the results are affected if 
we choose different values of A-p Q and -Efoid , which arc un- 
fortunately still uncertain even using the Suzaku data. We 
are aware that fixing these parameters is an oversimplica- 
tion because both A Fc and £f id could vary depending on 
the reflection parameter R. Using a fixed value of -Efdd = 
100 keV, both the reflection parameter and the power-law 
photon index are only slightly affected, where the best fit 
values provides 1.70 < T < 1.76 and 0.43 < R < 0.53, re- 
spectively for a range of 0.5 < ^Fe < 2.0. If we set Efo\d = 
50 keV, reflection parameter becomes slightly large (0.77 

< R < 0.97 for 0.5 < A Fc < 2.0), but statistically not sig- 
nificantly different. The power-law photon index is within 
a range of 1.70 < T < 1.76. Similarly, a different choice of 
E md = 150 keV provides 0.37 < R < 0.45 for 0.5 < A Fc 

< 2.0, with 1.72 < T < 1.76. Again these are relatively 
small effects and thus we assume A Fe = 1.0 and -Efoid = 
100 keV in the following analysis. 

4- 3. Iron Line Complex 

4.3.1. (1) Line Profile 

We next consider the X-ray spectrum of 3C 120 be- 
tween 3 and 10 keV, with the inclusion of the 5—7 keV 
data. The iron K line profile was then fitted in several 
steps. Firstly we fit the joint FI-XISs (XIS 0,2,3) and BI- 
XIS (XIS 1) hard X-ray spectra with the best-fit power- 
law plus reflection model (pexrav: PRV) as described 
in the previous section. The fit is very poor as shown 
in Figure 10 (a) and Table.3 (xl = 2500/1606 = 1.56; 
model- 1 in Table 3), mostly due to a prominent line pro- 
file near 6.2 keV in the observed frame (~6.4 keV in the 
rest frame). Then, adding a single Gaussian Fe K a line to 
the model (PRV+G) gives a greatly improved fit statistic 
(xl = 1807/1603 = 1.13; modcl-2). The line energy, con- 
verted into the rest frame, is 6.378to'oog keV, indicating 
neutral or low-ionization Fe. The measured width of the 
line is c b s = 193l'3 2 an d the equivalent width is EW 



= lOl^g eV. For completeness, we also added a small 
Compton shoulder to the iron K Q line, represented by a 
narrow Gaussian centered at 6.24 keV, with normalization 
fixed to 20% of the K a flux (Matt 2000), but there was no 
improvement to the fit at 90% confidence (no significant 
improvement in x 2 )- 

After the inclusion of the iron K line, there still re- 
mains a clear excess of photon counts in the Suzaku 
data between 5 and 7 keV (Figure 10 (b)). We fur- 
ther tried to improve the fit by assuming a rcdshiftcd 
disklinc rather than a simple Gaussian function (PRV 
+DL: modcl-3), or by adding another Gaussian to the 
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Fig. 10. Data/model ratio residuals of the 3 FI-XISs spec- 
trum of 3C 120 between 4—9 keV. upper (a): residual after 
fitting a power-law (r=1.73) with a reflection component (R 
= 0.47: model 1 in Table 3). middle (b): residuals after in- 
clusion of the narrow K Q core (model 2 in Table 3), bottom. 
(c): residuals after the addition of a broad diskline emission 
centered at 6.4 keV, as well as a 6.9 keV line (model 6 in Table 
3). 

model (PRV + 2G; modcl-4). Indeed the latter model 
(modcl-4) is suggested to fit the XMM- Newton data pretty 
well (Ogle et al. 2005), and also improves the x 2 statis- 
tic of the Suzaku data significantly (xl = 1757/1601 — 
1.10). However the inclusion of an additional Gaussian 
still does not represent very well the "red tail" below 6 
keV observed in the Suzaku data. Instead, inclusion of a 
broad diskline from an accretion disk (Fabian et al. 1989; 
see also Reeves et al. 2006 to fit the similar iron line pro- 
file of MCG-5-23-16) provides a better representation of 
the data. As a result, the fit statistic of our "best-fit" 
model (PRV+G+DL) was improved to xl = 1691/1599 
= 1.06 (model-5). This "PRV+G+DL" model (model- 
5), however, predicts a large inclination angle of i ~ 50 
deg. This seems problematic when taking into account 
the tight constraints on the jet inclination angle of i < 14 
deg from the superluminal motion, unless the jet is signif- 
icantly "warped" between the disk and the VLBI (i.e., pc 
scale) region. Also note that the centroid energy of the 
Fe line core is 6.368!q'qq5 keV, which is shifted redwards 
by ~ 30 eV compared to what is expected from a neutral 
iron K a line (6.400 keV). This shift is effectively larger 
than the uncertainties in the energy scale determined by 
the calibration sources (see § 2.1). 

We therefore considered whether we could fit the same 
line profile differently, to achieve a more face-on orienta- 
tion. Basically, the high inclination angle is driven from 
fitting the blue-wing of the line above 6.4 keV. If we as- 
sumed there is ionized emission from cither He-like (6.7 
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Table 3. Results of the spectral fits to the 3—10 kcV XIS 0—3 time-averaged, co-added spectrum with the different models used to 
describe the Fe-linc complex of 3C 120. 



Model 


Ei 


af 


E2 


°2 


Edl 


Tin 


i 


x 1 




[kcV] 


(EW) [eV] 


[kcV] 


(EW) [cV] 


[keV] 


[r g ] 


[deg] 


(d.o.f) 



1. PRV a 

2. PRV a +G b 

3. PRV+DL C 



a 070+O.OI3 
u -°'°-0.009 



1931;^ 
(101±iJ) 



fi onn+0.007 
o.o»u_ QQ 9 



40.7 



+6.'. 



18.2' 



4. PRV+2G 6.373±0.009 

5. PRV+G+DL 6.368±g:gg| 

6. PRV+2G+DL 6.398±0.01 



+ 15 
17 

+9 N 

+6 



158 
(90 
69 

(49±I) 
111+ 11 
(60±6) 



6.93±0.02 



6.93±0.02 



10 / 
(18) 



10 / 
(17) 



6.41 



+0.10 

-0.06 



19.4 



+5.1 



49.7 



+1.9 



6.40-f 



1.6 



+1.0 
-0.6 



6.5 



+3.8 
-3.2 



2500 
(1606) 

1807 
(1603) 

1853 
(1603) 

1757 
(1601) 

1691 
(1599) 

1681 
(1598) 



a pexrav model in XSPEC (Magdziarz & Zdziarski 1995). We assumed the direct power- law 
F3_iokeV = 3.0xl0 -11 erg cm -2 sec -1 (see Table 2). We fixed the reflection parameter R = 0. 
b A simple GAUSS function in XSPEC modified by redshift of 3C 120. 
c DISKLINE model assuming a power-law dependence of emissivity, /3 = —3.0. 
d An intrinsic width of iron line (excluding residual width in calibration <r ca i 
e 3-10 keV flux in units of 10" 12 erg cm" 2 s" 1 . 
f parameters fixed to these values. 



of T = 1.73 with 3-10 keV flux 
47 and Fe abundance Af c = 1.0. 



45±1 eV; see § 2.1). 



keV) or H-like (6.97 keV) Fe K a as was suggested in 
the Chandra HETGS and the XMM-Newton observations, 
that would reduce the inclination, as the blue-wing of the 
diskline would not need to extend much beyond 6.4 keV. 
While the fit with this "PRV+2G+DL" model (model-6 in 
Table 3) is only marginally better than before, the diskline 
now has a reasonable inclination of i = 6.5^|' \ deg with 
an EW of 32±5 eV, and fits the red-wing of the line profile 
well (while the ionized Fe emission models the blue-wing: 
see Figure 10(c)). It is also interesting to note that the 
centroid energy of the narrow Fe line core increases from 
6.368 keV to 6.398 keV (perfectly consistent with 6.400 
kcV within a ler error) in this face-on model (i < 10 deg). 

The measured width of narrow Fe line core (<7 bs ^120 
eV) is much broader than the residual width (er C ai) of a 
calibration line. Assuming a ca \ = 45 eV, the intrinsic 
width of Fe K Q line is cr int = U1±{1 eV (or EW of 60±6 
eV) . This line width is consistent with what was measured 
with XMM-Newton very recently by Ogle et al. (2005) and 
Ballantyne et al. (2004). Meanwhile, the additional broad 
disk line emission provides a better representation of the 
data if we fixed the center energy at 6.4 keV. Here we have 
assumed the outer radius of the diskline is r out = 1000 r g , 
and a typical steep emissivity of r -3 . The inner radius of 
the disk is constrained to be n n = 8.6^0 r g . A summary 
of the line parameters is given in Table 3. 
4-3.2. (2) Line Variability 

We next consider the iron line variability and its rela- 
tion to the variation of the baseline continuum component. 
Due to limited photon statistics, we traced the variabil- 
ity of the Fe K Q line core only, using a fixed line width 



and centroid energy. For this, a simple power-law plus 
Gaussian model is adequate to represent the Suzaku data 
in the 3—10 kcV band. Thus in our approximation, the 
only parameter which represents line variability is its nor- 
malization (i.e., integrated photon counts in iron line core) 
measured in different observation epochs. We are aware 
this is an oversimplified assumption given the complicated 
line profile and that variations with time are even possi- 
ble for the reflection component. Nevertheless, we believe 
this is the simplest way of searching for line variability on 
short timescalcs. The XIS data were divided into 4 equal 
time intervals (typically ~10 ksec) per each #1 — 4 obser- 
vations. Even assuming this simple PL + G model, all the 
datasets provide acceptable fits in the statistical sense of 

P(X 2 ) > io%. 

Figure 11 shows the spectral variability of the base- 
line power-law component and the iron K a core emission. 
Temporal evolution of 3—10 keV continuum flux, photon 
index, and Fe line core flux are plotted in the upper, mid- 
dle, and bottom panels, respectively. Although variability 
at some level may be present in the light curve of Fe line 
core flux, a constant fit provides a good explanation of the 
data, where x 2 /d.o.f = 17.3/15 (P(x 2 ) = 30%) and an av- 
erage iron line photon flux of F a = (4.41±0.19)x 10~ 5 ph 
cm" 2 s . A correlation between the underlying power- 
law continuum and Fe line fluxes is uncertain, but seems 
to be absent as shown in Figure 12. These results, how- 
ever, do not necessarily reject any variations of the iron 
line profile itself, if the intensity has not changed signifi- 
cantly during the observation. In fact, the iron line profile 
observed with Suzaku seems to have a somewhat differ- 
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Fig. 11. Spectral variability of baseline power-law compo- 
nent and iron K a core emission. 3—10 keV flux (upper), 
photon index (middle), and line flux (bottom) are separately 
shown. 4 XISs data are used for the fitting. 
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Fig. 12. Comparison of the iron K a core flux versus baseline 
continuum component. Variability is not clear for the iron line 
flux. 

ent shape than that observed with the XMM-Newton in 
2003. It seems that only the Suzaku data requires a strong 
red-tail, but instead the 6.9 keV emission seems to be a lit- 
tle more "spiky" in the XMM-Newton data (but we need 
more photon statistics to discriminate this further). A 
further detailed analysis of the line/absorption features is 
now on-going and will be presented in a forthcoming pa- 
per, especially by direct comparison of the XMM-Newton 
and the Suzaku spectra (Iwasawa ct al. in prep). 

4-4- Excess Soft X-ray Emission 

Finally we try to model the excess soft X-ray emission 
below 3 keV, as presented in Figure 7. We found both a 
steep power-law (r = 2.68±0.04) or thin thermal bremss 



emission (kT = 1.06±0.04 keV) can equally well fit the 
data with x 2 /dof ~ 1.1. Also the inclusion of an emission 
line feature near 0.9 keV improves the fit slightly, though 
it might an artifact of modeling the XIS data around Ne 

1 K edge (see Ogle et al. 2005 for XMM-Newton RGS soft 
X-ray spectrum). Interestingly the excess emission below 
3 keV seems a bit smaller in this Suzaku spectrum than 
that reported in the XMM-Newton observation (Figure 

2 of Ballantyne et al. 2005), for which the source was 
in a brighter state than any of the Suzaku observations 
(obs #1 — 4). Such excess emission was not necessary in 
the ASCA analysis (Grandi et al. 1997), but was found 
in both the ROSAT (Grandi et al. 1997) and BeppoSAX 
data (Zdziarski & Grandi 2001), as well as in the XMM- 
Newton data. 

This kind of soft excess emission has been observed in 
other BLRGs (Wozniak et al. 1998) and Sy-1 galaxies, and 
its nature is currently under debate. For 3C 120, it was 
first discovered in the Einstein SSS observations (Petre et 
al. 1984). Grandi etal. (1997) fit the 0.2-2.4 keV ROSAT 
data and obtained the best fit results with a steep power- 
law emission of T = 2.5—3.3, whereas Zdziarski & Grandi 
favor thin thermal emission of kT ~ 1 keV. Very recently, 
Ogle ct al (2005) suggested a steep power-law with T = 
2.7±0.1, which breaks at 0.6 keV to represent the XMM- 
Newton data below 1 keV. In contrast Ballantyne et al. 
(2004) assumed a thermal brcmsstrahlung component in 
order to reproduce the same X-ray data. Therefore the sit- 
uation is not still conclusive even after the XMM-Newton 
and the Suzaku observations. But if the soft excess is due 
to collisionally ionized plasma in an extended halo sur- 
rounding the nucleus as suggested in Zdziarski & Grandi 
2001, this should vary very little (but see the discussion 
in Ballantyne et al. 2004) . 

In fact, in some Seyfert galaxies, the soft excess emis- 
sion is consistent with being constant and is likely to 
be a combination of scattered power-law emission and 
photoionized gas associated with the AGN Narrow Line 
Region etc (e.g., Vaughan & Fabian (2004); Reeves et al. 
(2006)). Statistically wc cannot rule out there being some 
thcrmal/brcmss emission, but we will discuss below the 
power-law continuum emission which makes more sense 
in the context of the two components required to explain 
the spectral variability (see § 4.5). Moreover, with the 
observed variability being more strongly concentrated at 
lower energies, a more natural interpretation appears to 
be a steep power- law component. 

4-5. Multiband Spectral Variability 

Figure 13 shows the unfolded multiband spectra of 
3C 120 obtained with Suzaku in both high (obs #1) and 
low (obs #4) states. The best-fit model described in the 
previous sections is taken into account as various emission 
components in these figures; namely a direct power-law, 
Compton reflection, a narrow Fe line core, a 6.9 keV line, 
a broad disk line, and a steep power-law emission to rep- 
resent the soft excess emission below 3 keV. We first as- 
sumed the values determined in the previous sections as 
initial input parameters, and then re- fit the data again to 
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Fig. 13. Broad-band (0.6-50 keV: 3FI-XISs + PIN) spec- 
tra of 3C 120 observed with Suzaku in high (upper: obs #1) 
and low (lower: obs #4) states, with unfolded best-fit model 
components, red: total, blue: direct power-law component, 
green: reflection, purple: iron lines, red dashes; variable 
power-law component. The best fit parameters are summa- 
rized in Table 4. 

find a x 2 minimum under the constraint of fixed iron line 
parameters. This is because if we limit the analysis only 
above or below 3 keV as we have done in §4, this avoids 
having to include a scattered continuum component which 
could affect the fit. Therefore we thawed all the param- 
eters except for iron line parameters, to find a "new" \ 2 
minimum. The results of the fits are summarized in Table 
4 for various cases; (1) time-averaged spectra (obs #1 —4 
summed), (2) high flux state (HF; obs #1), and (3) low 
flux state (LF; obs #4), respectively. Direct comparison of 
the multiband spectra between HF and LF states (Figure 
14) clearly indicates energy dependence of spectral evolu- 
tion in different states of source activity. 

Thus it is interesting to consider what spectral compo- 
nent is primarily responsible for the spectral variability in 
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Fig. 14. Comparison of the best-fit broad-band (0.6—50 
keV) model of 3C 120 observed with Suzaku in high flux (HF; 
red: obs #1) and low flux (LF; blue: obs #4) states with the 
best-fit model functions. The best fit parameters for HF/LF 
states arc summarized in Table 4. 

3C 120. As can be seen in Table 4 and Figure 13, the 
underlying continuum emission (i.e., the sum of the direct 
power-law plus reflection) varies only a little. In fact, the 
spectral shapes, as measured by a photon index of T ~ 
1.75 in LF and HF states, are consistent within the sta- 
tistical error. The flux actually changes between LF and 
HF states, but only at the 10% level between 0.6—50 keV. 
In contrast, a factor of >5 variability has been observed 
for the steep power-law component where the 0.6—3 keV 
flux changes significantly from (9.3±1.8)x 10 -12 erg cm -2 
s _1 to (l.lj^'f)xlO -12 erg cm" 2 s" 1 . Actually, this steep 
power-law component is not visible when the source is in 
a low flux state, and only appears when the source is in 
a brighter state. This is consistent with the fact that the 
XMM-Newton spectrum shows a slightly larger amount of 
soft excess emission when the source was in a relatively 
bright flux state compared to any of the Suzaku observa- 
tions described in this paper. 

To further examine the origin of the spectral variabil- 
ity, the difference spectrum of HF minus LF spectra was 
extracted, using both the XIS and HXD/PIN data. The 
difference spectrum is plotted in Figure 15. This tech- 
nique clearly shows the variable component of the emis- 
sion from 3C 120 modified by absorption, with the con- 
stant components in the spectrum being subtracted. The 
resulting difference spectrum was fitted very well with a 
simple power-law of photon index T = 2.22±0.02 and fixed 
Ah = 1.23xl0 21 cm" 2 , with the fit statistic of x 2 /dof = 
198/192. Note that the inclusion of a thin thermal emis- 
sion component provides a rather poor fit as shown in 
Table 5. 

Although this spectral shape is a bit flatter than a steep 
power-law describing the soft excess emission (r ~ 2.7; 
see PL of Table 4), only small variations in the underly- 
ing continuum emission (PRV: r ~ 1.7) can easily account 
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Table 4. Results of spectral fits to the 0.6—50 keV XIS+PIN time-averaged, co-added spectrum of 3C 120. 



Component 


Parameter 


Average 


High (obs #1) 


Low (obs #4) 


wabs 


7V H a 


1.23' 




- 


PRV 


r 


1.74±0.02 


1.76±0.03 


1.78±0.03 




Sfoid [keV] 


100' 


— 


— 




R 


0.65±0.12 


0.79±0.22 


0.89±0.20 




i [deg] 


18.2' 


— 


— 




^0.6-3kcV 


26.2±0.5 


28.6±1.5 


26.5±1.0 




zp b 
^3-lOkoV 


28.7±0.2 


30.5±0.5 


27.6±0.3 




TP b 

10-50kcV 


63.3±0.8 


68.4±1.5 


61.8±1.3 


Gauss 1 


E [keV] 


6.40' 


- 


- 




a [eV] 


111' 


- 


- 




Nr C 
1 v line 


3.14' 


— 


— 


Gauss2 


E [keV] 


6.93' 


- 


- 




a [cV] 


10' 


- 


- 




Mi„c C 


0.69' 


- 


- 


Disklinc 


E [keV] 


6.40' 


- 


- 




Rm [r s ] 


8.6' 








i [deg] 


6.5' 








1 v line 


1.57' 






PL 


r 


2.66±0.04 


2.60±0.08 


2.65' 




-^0.6-3keV 


4.8±0.6 


9.3±1.8 


1 1+ 13 
1 ' L -i.i 




TP b 
^3-10kcV 


1.4±0.2 


2.9±0.6 


3 +0 - 4 


Gauss3 


E [keV] 


0.87' 








a [eV] 


10' 








iV linc c 


5.71' 






xVd.o.f 




2978/2768 


3018/2884 


2853/2729 



a Galactic absorption column density in units of 10 21 cm -2 . 

6 Flux in units of 10" 12 erg cm" 2 s — . 

c Line normalization in units of 10~ 5 photons cm -2 s" 1 . 

Assumed to be same as the left, i.e., those assumed to fit an "average" spectrum. 
^ Parameters fixed to these values. 



Table 5. Results of spectral fits to the 0.3—50 keV difference spectrum between obs #1 and obs #4. 



Model 


kT [keV] 


brm 


n 




1 2 


^PL2 


X 2 /d.o.f 


wabs + Brcmss 


2.53±0.07 


13.0±0.5 










453/192 


wabs + PL 






2.22±0.02 


14.8±0.2 






198/192 


wabs + 2PL 






2.65' 


9.4±0.3 


1.75' 


6.3±0.3 


198/192 



a A thin-thermal or power-law function modified by Galactic absorption of Nh = 1.23xl0 23 cm 2 . 
b 0.6-10 keV flux in units of 10" 12 erg cm" 2 s" 1 
c Differential spectral photon index. 
* Parameters fixed to these values. 
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Fig. 15. The difference spectrum of 3C 120 between high 
(obs #1) and low (obs #4) flux states in the Suzaku observa- 
tion (3 FI-XISs + PIN). Although the difference spectrum is 
well represented by an absorbed power-law of T = 2.22±0.02 
in the XIS energy band, some residuals remains above 10 keV. 
The overall spectrum (including HXD/PIN) is better repre- 
sented by a double power-law function, where the low energy 
emission (below 3 keV) is dominated by a steep power-law of 
r = 2.65 (red dashes) and the high energy emission (above 
10 keV) is modified by a flat power-law of V = 1.75 (blue 
dashes). 

for this apparent discrepancy. In fact, the same difference 
spectrum can be equally well fitted with double power-law 
functions, where the low energy emission is dominated by 
a steep power-law (ri = 2.65, which mimics the soft "vari- 
able" power-law) and the high energy emission is repre- 
sented by a flat power-law (T2 = 1.75, which mimics the 
hard "direct" power-law component). These model com- 
ponents are presented in Figure 15 as red and blue dashes, 
respectively. 

We also note that a normalization of the HXD/PIN 
difference spectrum seems a bit larger than an extrapola- 
tion of the XIS spectrum above 10 keV. This may imply 
both intrinsic and reprocessed component change together 
in 3C 120, though it was not the case for MCG-6-30-15, 
where the reflected component is constant and only the 
intrinsic power-law varies (Minuitti et al. 2006: see also 
§5.3). But we suspect a part of the reason could be in- 
strumental: As we see in §2 and in the Appendix, current 
background models of the HXD/PIN would produce ~5% 
systematics in the background subtraction, which corre- 
sponds to ^20% uncertainties in the 3C 120 flux above 
10 keV. This could be more enhanced in the differnce 
spectrum as shown in Figure 15. We will discuss more 
about this in a forthcoming paper using update back- 
ground models of the HXD/PIN. 



5. Discussion 

5.1. Suzaku View of 3C 120 

In previous sections, we presented temporal and spec- 
tral analysis of a Suzaku observation of 3C 120 in February 
and March 2006. The great advantage of Suzaku is that 
both the XIS and HXD/PIN have excellent sensitivity, so 
that we can resolve the spectral evolution of this source 
over two decades in frequency even within a short (40 ksec) 
exposure. During the Suzaku observations, 3C 120 was in 
a relatively low flux state of i^-iokcV = 3.9xl0~ n erg 
cm~ 2 s _1 on average, though the flux decreased gradu- 
ally from 4.4x 10~ n erg cm~ 2 s _1 to 3.6x 10~ n erg cm -2 
s _1 between obs #1 and obs #4. Historically, the maxi- 
mum flux observed by Suzaku (that obtained in obs #1) 
is almost consistent with what was observed by ASCA 
in 1994, by Beppo-SAX in 1997 and by XMM-Newton 
in 2003 (iViokoV ~ 4.6xl0" n erg cm' 2 s" 1 ; Reynolds 
1997; Zdziarski & Grandi 2001; Ballantyne et al. 2004) 
but 40 % lower than that reported by RXTE in 1998 
(Eracleous et al. 2000). This range of variability is nat- 
ural considering the long term variability of this source 
over year-long timescales in X-rays (Marscher et al. 2002; 
Markowitz & Edelson 2004). 

Using the high sensitivity, broad-band instruments on- 
board Suzaku, we have confirmed various important char- 
acteristics of 3C 120 ; (1) the presence of a relatively nar- 
row (cj = mtU cV or EW of 60±6 eV) Fe line core 
centered at 6.4 keV (see Ballantyne et al. 2004; Ogle et al. 
2005 for XMM-Newton analysis). In addition, Suzaku has 
unambiguously confirmed the broad iron line component 
which may be possibly related with an extremely broad 
line feature observed with ASCA (Grandi et al. 1997; 
Reynolds 1997; Sambruna et al. 1999), (2) a relatively 
weak reflection component (i? = 0.65±0.12: see Table 4), 
determined with much better accuracy that observed with 
Beppo-SAX (R = 0.56t^ u ; Zdziarski & Grandi 2001) 
and RXTE (R = 0A±%i; Eracleous et al. 2000; sec also 
Ballantyne et al. 2004 for updated results) under the same 
constrained parameters, (3) the presence of a soft excess 
emission component below 3 keV, as firstly reported by 
Einstein (Turner et al. 1991) and by ROSAT (Grandi et 
al. 1997). Finally, we confirmed the spectral variability 
of 3C 120 whereby the spectra becomes softer when it be- 
comes brighter, as suggested by Maraschi et al. (1991) and 
Zdziarski & Grandi (2001). 

The important question possibly raised by readers is 
"what is completely NEW for Suzaku!" . Firstly, we 
discovered, for the first time, the presence of a broad 
red-tail (asymmetric red wing) in the Fe K Q line of any 
radio galaxy. Interestingly, the line profile of 3C 120 
is quite similar to those found in some Seyfert galaxies 
(e.g., MCG-5-23-16; Reeves et al. 2006), which provides 
important challenges to the unification models of radio- 
loud/quiet AGNs. We also detected a 6.9 keV line, but 
perhaps not as clearly as first noticed in the XMM-Newton 
data in 2003 (Ballantyne et al. 2004). This may suggest 
significant changes in the Fe line profile on the year-long 
timcscalc, but further deep/long monitoring of 3C 120 is 
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necessary for this kind of study. Also, we have shown 
that Suzaku can provide a new diagnosis in accurate mea- 
surement of the line center energy (§4.3), which was not 
possible before. 

Second, we found the excess variance (i.e., variability 
amplitude) is generally larger at lower energies. From de- 
tailed multiband spectral studies, we conclude the primary 
"variable" component is a steep (r = 2.6—2.7) power-law, 
which, at the same time, accounts for the soft excess emis- 
sion below 3 keV. In fact, the difference spectrum is well 
represented by a steep power-law (ri = 2.65) modified by 
a flatter second component (I^ = 1.75) at higher ener- 
gies. Apparently, the flat component mimics the "direct" 
power-law emission from the 3C 120 nucleus, and varies 
only a little (10 %) during our observations. In contrast 
the flux of the steep power-law emission may have changed 
by a factor of > 5 in the transition between the HF and LF 
states. Below we will discuss these new findings provided 
by Suzaku and their interpretation in detail. 

5.2. The Nature of the Iron K Line Complex 

Suzaku has successfully resolved the iron K line complex 
of 3C 120 and also was first to verify the broad compo- 
nent's asymmetry. We showed that the iron line complex 
is composed of (1) a relatively narrow, neutral iron K line 
core, (2) broad iron line emission possibly emitted from 
the accretion disk (see § 4.3), and (3) an ionized ~6.9 
keV line. Interestingly, any of the models assuming a sim- 
ple Gaussian profile (modcl-2, 4, 5 in Table 3) and/or an 
"edge-on" disk (model-5) predicts that the centroid en- 
ergy of the narrow line core is slightly shifted redwards 
by ~ 30 cV. Meanwhile, the energy of the narrow core in- 
creases from 6.368 keV to 6.398 keV in the face-on model 
(model-6), which is then formally consistent with neutral 
iron at 6.40 keV. We suspect this is because the face-on 
diskline contributes towards some of the flux of the nar- 
row 6.4 keV line core. For instance, the observed 30 eV 
shift of the line core implies a typical radius of ~200 r g if 
the red-shift is purely gravitational. 

In contrast, the intrinsic line width of the narrow Fe 
line core was measured to be a ~ 110 eV. Assuming this 
line width is simply caused by Doppler broadening, this 
corresponds to a FWHM velocity of 10 4 km s -1 , indicating 
a possible origin in the broad line region (BLR). In the 
optical, the FWHM velocity of the BLR, as measured by 
Peterson et al. (2004) using the Hg is smaller, 2200±200 
km s _1 . Alternatively, Ogle et al. (2005) provides a very 
high S /N optical spectrum which shows a broader (3> 2000 
km s _1 ) component associated with the BLR lines, which 
would be then more consistent with the width of the Fe 
line core if that also originates from the BLR. Thus we 
suggest most of the flux of the narrow 6.4 keV line core 
could be from the BLR, but that the outer disk (~200 r g ) 
also makes some contribution. 

The residuals present after subtracting the iron line core 
(Figure 10 (b)) are poorly modeled either by adding a sim- 
ple broad Gaussian or by adding the Compton shoulder of 
the K a line, as we have discussed in §4.3. One may also 
suggest whether the presence of a warm absorber may 
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affect continuum curvature to make an apparent red-tail 
below the iron K line. In this context, Ogle et al. (2005) 
present a very high S/N RGS spectrum, and there do not 
appear to be any significant lines or edges due to a warm 
absorber. So the warm absorber must be very weak or 
even absent in 3C 120 (but not surprising if we are view- 
ing face-on). This also means that it is rather unlikely 
that any additional absorption can effect the Fe K line 
modeling, especially the red-wing (see also §4.5). Rather, 
a significant red-tail below 6.4 keV (in the rest frame of 
source) favors diskline emission from the inner accretion 
disk of n n ~ 8.6^Qg r g . Indeed, adding the diskline emis- 
sion provides a much better fit statistic than any other 
model as summarized in Table 3. If this broad line re- 
ally originates from the inner accretion disk, it provides 
important clues to jet formation in the accretion disk. 

For example, Reeves et al. (2006) have discovered a sim- 
ilar broad iron line in the Suzaku/XMM-Newton spectra 
of MCG-5-23-16, which is thought to originate from inner 
accretion disk (rj n ~ 20 r g ). These observations may im- 
ply that both the radio-loud 3C 120 and the radio-quiet 
MCG-5-23-16 have similar accretion disk structure, in 
contrast to suggestions that the optically-thick accretion 
disk is truncated in 3C 120 to a hot, optically thin flow at 
a distance of r ln ~ 100 r g (Eraclcous et al. 2000; Zdziarki 
& Grandi 2001; Ballantyne et al. 2004). We also note 
that an EW of the diskline is 32±5 eV, whereas that for 
the narrow Fe line core is 60±6 eV. Given the prediction 
of George and Fabian (1991) that a reflection parameter 
R should be equal to EW(eV)/150(eV), the iron line EW 
is expected to be 90 eV for the level of Compton reflec- 
tion R ~ 0.6. Interestingly, this is consistent with the sum 
of the EWs of the narrow and diskline components, and 
therefore the Compoton hump could possibly be associ- 
ated with both lines. Observations of the iron line profiles 
in various other broad line radio galaxies are important for 
the systematic comparison between Seyferts and BLRGs. 
We are planning to submit further deep observations of 
other BLRGs in the next Suzaku observation program to 
test this. 

Finally, to further test the robustness of the iron K 
diskline, we tried an alternative partial covering model, to 
see if that can reproduce the broad residuals present in the 
iron K band. Instead of the broad line, the neutral par- 
tial covering XSPEC model pcfabs was included, and the 
XIS+HXD data re-fitted over the range from 0.6—50 keV. 
The model assumes an absorbed double power-law con- 
tinuum, together with a single Gaussian line to model the 
iron Ka line core, and ionized iron K emission as described 
previously, while Compton reflection is also included in 
the model with a cut-off energy of 100 keV. The par- 
tial covering model results in a fit statistic of 2999/2769 
and leaves significant residuals around 6 keV in the iron 
K band. The fit statistic is statistically worse than the 
diskline fit (A^ 2 = 21 for 1 extra parameter), while the 
diskline has fit statistic of 2978/2768 (see §4.5), and the 
equivalent width of the diskline- is 45 eV. Without includ- 
ing cither a diskline or a partial coverer, the fit statistic is 
3018/2771. For a 50% covering fraction, the upper-limit 
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Fig. 16. Comparison of the multiband spectrum of a ra- 
dio-loud BLRG 3C 120 (upper) and a radio-quiet Sy-1 galaxy 
MCG-6-30-15 (lower). In both figures, X-ray data comes 
from the most recent Suzaku observations (3C 120: this work. 
MCG-6-30-15: Minuitti et al. 2006) and the other data comes 
from the NED data base. Red lines show the best-fit X-ray 
spectral model during the high state, whereas blue lines shows 
that for the low state. Variable power-law component is 
shown as a thick dashed line for 3C 120, whereas the differ- 
ence spectrum is shown for MCG-6-30-15. The thin dotted 
line corresponds to an example fit of the one-zone homoge- 
neous SSC model as described in the text. 

on the column density is then < 3 x 10 21 cm~ 2 . Therefore 
the detection of the broad iron line appears to be robust 
in 3C120. 

5.3. "Hidden' Jet Emission? 

In § 4.5, we considered the origin of the spectral vari- 
ability in 3C 120 by comparing the multiband spectra 
obtained in low flux (LF) and high flux (HF) states. 
Interestingly, several authors have examined similar spec- 
tral evolution of radio-quiet Seyfert galaxies. For exam- 
ple, Miniutti et al. (2006) compared multiband spectra of 
the Sy-1 MCG-6-30-15 obtained in the LF and HF states. 
They found that difference spectrum shows a steep power- 



law of F = 2.2, which is consistent with the direct power- 
law emission observed in both LF and HF states (see also 
Figure 16). This leads to an idea that the broad band 
spectral variability of MCG-6-30-15 is decomposed into 
two components; a highly variable power-law (direct nu- 
cleus emission) and a constant reflection component plus 
iron line emission. Similarly, Reeves et al. (2006) pre- 
sented a difference spectrum of a Sy-2 galaxy, MCG-5- 
23-16, which was fitted extremely well by a simple ab- 
sorbed power-law of V = 1.9. Again, this spectral shape 
is exactly consistent with the direct power-law component, 
which they observed in LF and HF states. Furthermore, 
there are no residuals present in the iron K band or any 
excess counts in the HXD/PIN difference spectrum above 
10 keV. 

In contrast to these findings for Seyfert galaxies, an in- 
teresting discovery made by Suzaku is that the variable 
component in 3C 120 is much steeper (T ~ 2.7) than the 
power-law emission reported in literature (1.6 < Y < 1.8: 
see also Marscher et al. 2002). One interesting idea to 
account for this steep, variable emission is the beamed 
radiation from the jet (not on kpc scales as discussed in 
Harris et al. 2004 and §1, but the unresolved base of the 
jet on sub-pc scales), though this component contributes 
only ~ 20% at most, of the Sy-1 like X-ray emission in 
3C 120 (i.e., emitted from the disk and corona). In fact, 
3C 120 possesses a supcrluminal radio jet with an incli- 
nation angle < 14 deg (Eracleous & Halpern 1988). This 
low inclination angle implies that 3C 120 may have some 
"blazar-like" characteristics, such as rapid X-ray variabil- 
ity or a non-thermal spectrum extending to the 7-ray en- 
ergy band. For example, if we assume a jet bulk Lorentz 
factor of Tblk = 10 and a jet inclination angle of 14 deg, 
we expect the observed flux is mildly "boosted" towards 
the observer with a Dopplcr beaming factor 6 ~ 3. Note 
also that the X-ray spectra of blazars are generally rep- 
resented by a power-law function, but their X-ray photon 
indices ranges widely (1.5 < T < 3.0; Kubo ctal. 1998) 
among various sub-classes (e.g., HBL, LBL, and FSRQ; 
e.g., Fossati et al. 1998). Interestingly, similarly variable, 
steep power-law emission has been found for a distant 
(z = 0.94) quasar PG1407+265 which is classified as a 
radio-quiet object, but shows strong evidence of a radio 
jet with a highly relativistic speed in the VLBA observa- 
tions (Gallo 2006). 

Figure 16 shows the multi-band spectrum of 3C 120 
(upper), compared to that of the radio-quiet Sy-1 galaxy 
MCG-6-30-15 (lower) between the radio and X-ray en- 
ergy bands. Note the significant difference in the radio 
band (as inferred by a definition of radio-quiet/loud ob- 
jects), but spectral similarity in the optical band. In both 
figures, the X-ray data comes from the recent Suzaku ob- 
servations. Red/blue lines show the best-fit X-ray spectral 
model during HF/LF states. The difference spectrum is 
shown as a thick dashed line. The origin of the variable 
component (steep power-law) in 3C 120 is still uncertain, 
but here we assume the non-thermal jet dominates this 
power-law emission. The thin dotted line shows an exam- 
ple fit with a one-zone homogeneous SSC (Synchrotron 



Suzaku Observations of 3C 120 



17 



self-Compton) model described in Kataoka et al. (1999), 
where the soft X-ray excess corresponds to the highest 
energy end of the synchrotron emission. 

Also above 10 keV, the inverse Compton component is 
possibly present even in the hard X-ray band, but the 
isotropic Seyfert-like emission overwelms that blazar-like 
emission, where the relativistic boost of the jet is smaller 
than in blazars, because of a lower S. Note that this model 
for the jet-like emission - with the SED peaking at ~ 3 
x 10 15 — 10 16 Hz - implies HBL-like classification of the 
jetted emission. This in principle fits well with the FR- 
I type classification of the object discussed in Section 1, 
where HBL blazars are FR-I radio galaxies viewed close 
to the direction of the jet. However, there is an additional 
complication: the FR-I / HBL nature of this object as 
inferred from the radio observations and the soft X-ray 
excess is in conflict with the presence of broad emission 
lines, present in the optical spectrum, as those are gener- 
ally not observed in FR-I sources. 

Interestingly, the jet parameters derived here are consis- 
tent with typical values for blazars, except for a moderate 
beaming factor of 6 = 3 (viz 5 < S < 30 for blazars; e.g., 
Kubo et al. 1998): magnetic field B = 0.3 G, region size R 
— ct var 6 ~ 1.0xl0 17 cm, under an assumption of equipar- 
tition between electron and field energy densities (u c = 
ub)- We assumed a broken power-law form of the electron 
distribution N(j) cx 7~ s (l+7/7br) _1 cxp(-7/7 max ), 
where 7 is the electron energy (in units of m e c 2 ), 7t> r k 
= 10 3 , 7max = 1.3xl0 5 , and s = 1.8. Surprisingly, this 
simple jet model reproduces the general trend of the spec- 
tral energy distribution quite well, though some discrep- 
ancies are seen in the radio band. Such discrepancies in 
the radio are found in many studies and it is believed to 
be a consequence of the radio emission originating from a 
much larger region than the X-rays (e.g., Kataoka et al. 
1999). Also note that the prediction of gamma-ray flux 
via inverse Comptonization is less than 10 -10 erg cm~ 2 
s , qualitatively consistent with the non-detection by 
EGRET and COMPTEL onboard CGRO. 

Based on these facts alone, it is still premature to con- 
clude that the steep variable X-ray emission is actually 
originating from the jet in 3C 120. Maraschi et al. (1991) 
and Zdziarski & Grandi (2001) claimed that an X-ray jet 
component could, at the very least, dilute the Seyfert-like 
spectrum of 3C 120, and may account for the observed 
weak reflection features. This idea is straightforward, but 
may be oversimplified by following reasons: (1) the Suzaku 
observations have now revealed that the variable emission 
can explain at most 20% of the Sy-1 emission even be- 
low 3 keV. (2) The putative jet spectrum is much steeper 
than the direct power-law (r ~ 1.7). Such a low ampli- 
tude, steep power-law jet component cannot dilute the Sy- 
1 emission sufficiently, especially above 5 keV (where the 
iron lines and reflection become important, unless there is 
additional, weak contribution from a hard spectrum due 
to the inverse Compton component. 4 (3) No evidence 

4 In this context, we note that an inverse Compton model in 
Figure 16 is dependent on input parameters as fic oc u~ 2 , where 



for "blazar-like" variability has been found above 2 keV 
(Giozzi ct al. 2003; Marshall ct al. 2003). Therefore, it 
seems that the weak iron line and the reflection hump ob- 
served in BLRGs (Eracleous et al. 2000) might be more 
intrinsic, and provides important clues to the origin of the 
jet-like emission. 

Nevertheless, we showed that the variable steep X-ray 
emission is a key to understanding the spectral evolution 
of 3C 120. In particular, it seems that the steep power-law 
component is more significant when the source is brighter, 
and can well explain the significant soft excess repeatedly 
observed in the literature for this AGN. Interestingly, the 
soft excess emission is hardly visible when the source be- 
comes fainter, which may indicate that the jet component 
has completely disappeared and/or is hidden behind much 
stronger Sy-1 emission. Unfortunately, the flux changes 
of 3C 120 during the Suzaku observation were relatively 
small, so we cannot conclude what fraction of X-ray flux 
is actually explained by the jet. Future deep Suzaku ob- 
servations, as well as continuing VLBI monitoring coinci- 
dent with X-ray monitoring (as the campaigns reported 
by Marscher et al. 2002), sensitive measurements with 
GLAST, in quite different states of source activity (i.e., 
observations at historically high and low states) will be 
crucial to understanding the nature of emission proper- 
ties in 3C 120 (see recent paper by Grandi and Palumbo 
(2006) for the detectability with GLAST at MeV-GcV 
energy band). 

6. Conclusion 

We have presented a detailed analysis of results from 
the broad line radio galaxy 3C 120 observed with Suzaku 
from February and March 2006. Thanks to the excellent 
sensitivity of both the XIS and HXD/PIN detectors 
onboard Suzaku, we obtained multi-band data with 
unprecedented accuracy between 0.6 and 50 keV. Our 
major findings are as follows: 

(1) An overall spectral shape, including weak reflection 
(R ~ 0.6) and a soft excess which was consistent with 
previous findings, and was determined with unprece- 
dented accuracy. 

(2) We confirmed the presence of a narrow iron K Q core 
with a width of a = 111±H eV (and an EW = 60±6 eV), 
and a 6.93±0.02 keV emission line possibly emitted from 
H-likc iron. 

(3) After subtracting these narrow line components, 
significant residuals remain redwards of the narrow iron 
line core, well below 6.4 keV. 

(4) The difference spectrum between LF and HF clearly 
shows a power-law of V = 2.2, which is naturally in- 
terpreted as a combination of a highly variable steep 
power-law (r = 2.7; jet-like) plus a moderately variable 
direct power-law emission (r = 1.7; Sy-1 like). No 
significant variability was found in either the Fe K line 

flQ is an observed inverse Compton flux. Therefore this could 
dilute the Sy-1 emission only if the magnetic field strength is 
about a factor of five smaller than equipartition values, which 
is often not the case for blazars (see, e.g., Kubo et al. 1998). 
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emission or reflection component. 

We argue that the narrow iron K line is primarily emit- 
ted from the broad line region corresponding to a FWHM 
velocity of 10 4 km/s, but the emission from the outer disk 
(~ 200 r g ) may also play a part. Meanwhile, the broad 
redshifted iron line can be interpreted as emission from 
the inner, face-on accretion disk (r- m ~ 10r ff and i < 10 
deg). We considered the origin of a highly variable steep 
power-law component in the context of a possible rela- 
tion to non-thermal jet emission. Although we cannot 
conclude what fraction of X-ray emission in 3C 120 is ac- 
tually explained by the jet, it must be less than 20 % of 
the "Sy-1 like" emission even below 3 keV. This clear dif- 
ference of spectral evolution between Sy-ls and 3C 120 
provides an important motivation for further future deep 
observations of BLRGs. 

Finally, we thank an anonymous referee for his/her con- 
structing comments which helped clarify many of the is- 
sues presented in this paper. We also thank Dr. Tahir 
Yaqoob for his helpful comments and discussion on the 
Suzaku data analysis. 

Appendix 1. Effect of Background Variations on 
the HXD/PIN Light Curve 

Here we provide a follow-up discussion concern- 
ing the robustness of background subtraction of the 
HXD/PIN using the most recent response and back- 
ground model as of November 2006 (vl.2 data and 
ae_hxd_pinxinom_20060814.rsp). Detailed studies of back- 
ground systematics are still under investigation by the 
HXD instrumental team, but careful comparison of the 
HXD/PIN light curves (see Figure 3) with background 
variations may provide important hints for future model- 
ing of relatively faint sources. For this aim, 3C 120 is a 
good target as it is bright enough to be firmly detected 
by the HXD/PIN, but much fainter than the background 
photon statistic. As we have seen in Figure 2, the net 
intensity of 3C 120 is expected to be ~30 % of the non 
X-ray background plus CXB above 12 keV. 

Figure 17 compares temporal variations of total 
HXD/PIN counts (blue: including source and back- 
ground), non X-ray background model A (green; see §2.2), 
and resultant estimate of source counts from 3C 120 (red). 
The HXD/PIN count rate (blue) is highly variable from 
0.5 counts s _1 to 0.9 counts s _1 , and this trend is generally 
well reproduced in the background model (green). The 
dead time is corrected for each of the time bins, because 
it varies with the event rate (dominated by background 
photons) of the HXD/PIN detector and ranges from 3.5% 
to 8.5% during 3C 120 #1 — 4 observations. As a result, 
small residuals (of 5—10 % background level) remains after 
subtracting the background photons, especially when the 
total HXD/PIN counts exceeds ~ 0.7 counts s _1 . Since 
these time regions correspond to orbits where the COR is 
low and/or passing through the SAA, a long-decay com- 
ponent of the background may have not been modeled 
perfectly in the current data. 
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Fig. 17. Comparison of temporal variation of total 
HXD/PIN counts (blue: including source and background), 
non X-ray background model A (green; see §2.2), and 
resultant estimate of source counts from 3C 120 (red) for 
3C 120 #1-4. 
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Similarly, the relative lack of variability in the XIS light 
curve suggests apparent "rapid" changes in HXD/PIN 
flux Figure 17 could be artifacts of the background sub- 
traction. Nevertheless, the general trend of the HXD/PIN 
flux variations traces well those of the XIS data be- 
low 10 keV, as we have seen in Figure 5. The average 
net source count rate of the HXD/PIN is 0.195±0.005 
cts s" 1 , 0.177±0.005 cts s"\ 0.180 ±0.004 cts s" 1 , and 
0.165±0.004 cts s _1 respectively, and hence the maximum 
is in obs #1 and the minimum in $4, as per the XIS. A 
constant fit to overall HXD/PIN light curves (obs #1 — 4) 
provides x 2 /d.o.f = 154/80. Even after removing the time 
regions where the background is high (HXD/PIN counts 
exceeds 0.7 counts s _1 ), we obtain x 2 /d.o.f = 103/71 
(P{x 2 ) = 0-7 %). Therefore we think the long-term varia- 
tions of 3C 120 are robust in the 12—40 keV energy band. 
More detailed discussion using revised background models 
will be discussed elsewhere. 
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